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ABSTRACT 

We report the discovery of one of the most powerful sunquakes detected to date, 
produced by an X1.2-class solar flare in active region 10720 on 2005 January 15. We 
used helioseismic holography to image the source of seismic waves emitted into the 
solar interior from the site of the flare. Acoustic egression power maps at 3 and 6 mHz 
with a 2 mHz bandpass reveal a compact acoustic source strongly correlated with 
impulsive hard X-ray and visible-continuum emission along the penumbral neutral 
line separating the two major opposing umbrae in the (5-configuration sunspot that 
predominates AR10720. At 6 mHz the seismic source has two components, an intense, 
compact kernel located on the penumbral neutral line of the (5-configuration sunspot 
that predominates AR10720, and a significantly more diffuse signature distributed 
along the neutral line up to ~15 Mm east and ^--^30 Mm west of the kernel. The 
acoustic emission signatures were directly aligned with both hard X-ray and visible 
continuum emission that emanated during the flare. The visible continuum emission is 
estimated at 2.0 x 10^^ J, approximately 500 times the seismic emission of ~ 4 x 10^" J. 
The flare of 2005 January 15 exhibits the same close spatial alignment between the 
sources of the seismic emission and impulsive visible continuum emission as previous 
flares, reinforcing the hypothesis that the acoustic emission may be driven by heating 
of the low photosphere. However, it is a major exception in that there was no signature 
to indicate the inclusion of protons in the particle beams thought to supply the energy 
radiated by the flare. The continued strong coincidence between the sources of seismic 
emission and impulsive visible continuum emission in the case of a proton-deficient 
white-light fiare lends substantial support to the "back - warming" hypothesis, that 
the low photosphere is significantly heated by intense Balmer and Paschen continuum- 
edge radiation from the overlying chromosphere in white-light fiares. 
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1 INTRODUCTION 

Although most large solar flares appear to be acoustically 
inactive, certain energetic flares radiate intense seismic tran- 
sients into the solar interior during the impulsive phase. 
These wave packets radiate thousands of kilometres from the 
flaring region into the solar interior, but most of this energy 
is refracted back to the solar surface within approximately 
50 Mm of the source and within an hour of the beginning 
of the flare. The surface manifestation is a wave-packet of 
ripples accelerating outward from the general source region 
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that is sometimes obvio u s in r aw helioseismic observations. 
iKosovichev fc Zharkoval (Il998l ) discovered the first known 
instance of seismic emission, from the X2-class fiare of 1996 
July 9 in AR7978, identifying the phenomenon by the name 
"sunquake." 

For a long time these events were thought to be an ex- 
tremely rare phenomenon. However, with the advancement 
of loca l helioseismic techniques such as helioseismic holog- 
raphy (iLindsev and Braun|[200ol ). we have now detected nu- 
merous seismic sources of varyi ng intensity pr oduced by 
X- and high M-class flares (.Donea. Braun fc Lindsev, .19991 : 
iDonea fc Lindsevll2005l : lOonea et al.ll2006l ). 

A subsequent extensive survey of X-class solar fiares 
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(Besliu-Ionescu et al., in preparation) led to the discovery 
of more than a dozen seismic emission signatures from flares. 
Almost all of these have occurred in complex active regions. 
In this paper we report on the discovery of one of the most 
powerful flare seismic transients detected to date and com- 
pare the acoustic signatures of this sunquake with other sup- 
porting observations. 

AR10720 was a complex active region that appeared on 
the solar disk on 2005 January 11 and soon became one of 
the largest and most active sunspot regions of the current 
solar cycle. In the period January 15 - 20, AR10720 pro- 
duced 5 X-class solar flares, including an X7.1 on January 
20, which produced an intense solar proton storm. However, 
helioseismic observations sufficient to show seismic emission 
were acquired only for the XI. 2 flare of January 15. This flare 
was observed by numerous space and ground-based solar ob- 
servatories, including the Michelson Doppler Imager (MDI) 
instrument on board the Solar and Heliospheric Observa- 
tory (SOHO), the Reuven Ramaty High-Energy Solar Spec- 
troscopic Imager (RHESSI), the Geostationary Operational 
Environmental Satellite (GOES), the Transition Region and 
Coronal Explorer (TRACE), and the earth-based Global 
Oscillations Network Group (GONG). AR10720 itself was 
observed by the Imaging Vector Magnetograph (IVM) at 
the Mees Solar Observatory in the general time frame of the 
15 January 2005 flare. 

In Section 2 we present a brief description of the helio- 
seismic observations we analyzed and in Section 3 we review 
the technique of computational seismic holography. Section 

4 consists of our results and analysis, and flnally in Section 

5 we present a discussion and summary. 



2 THE HELIOSEISMIC OBSERVATIONS 

The MDI data consist of full-disk Doppler images in the pho- 
tospheric line Ni I 6768 A, obtained at a cadence of 1 minute, 
in addition to approximately hourly continuum intensity im- 
ages and line-of-sight magnetograms. The MDI d ata sets are 
described in more detail by Schorrc r et al.l l| 19951 . For the 
flare of 2005 January 15, we analyzed a dataset with a period 
of four hours encompassing the flare. For the purpose of our 
analysis, the MDI images we obtained (Dopplergrams, mag- 
netograms and intensity continuum) were remapped onto 
a Postel projection (jDoForcst 2004 1 that tracks solar ro- 
tation, with the region of interest fixed at the centre of 
the projection. The nominal pixel separation of the projec- 
tion was 0.002 solar radii (1.4 Mm) with a 256 x 256 pixels 
fleld of view, thus encompassing a region of approximately 
360 X 360 Mm^ on the solar surface. 



3 HELIOSEISMIC HOLOGRAPHY 

We briefly review the adaptation of computational seismic 
holography for applications in flare seismology. In general, 
helioseismic holography is the phase-coherent reconstruction 
of acoustic waves observed at the solar surface into the so- 
lax interior to render stigmatic images of subsurface sources 
that have given rise to the surface disturbance. Because the 
solar interior refracts down-going waves back to the sur- 
face, helioseismic holography can likewise use observations 



in one surface region, the pupil, to image another, the focus, 
a considerable distance from the pupil. We call this "seis- 
mic holography from the subjacent vantage" (see flg. 4 of 
iLindsev and BraunI (|2000l ')'). The subjacent vantage renders 
the photosphere as viewed by an acoustic observer directly 
beneath it. 

In general the acoustic reconstruction can be done ei- 
ther forward or backward in time. When it is backward in 
time, we call the extrapolated fleld the "acoustic egression." 
In the case of subjacent vantage holography, this represents 
waves emanating from the surface focus downward into the 
solar interior. When the surface acoustic field at any point 
r' in the pupil is expressed as a complex amplitude tp for any 
given frequency w, the acoustic egression can be expressed 
as 
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In this formalism, G+(r, r', w) is a Green's function 
that expresses the disturbance at the focus, r, due to a mea- 
sured point source at surface point r' from which the acous- 
tic wave is supposed to propagate backwards in time to the 
focus. 

The relation between the complex amplitude, ip{r, u), 
of frequency appearing in equation ([1} and the real acoustic 
field, ip{r, t), representing the surface acoustic field in the 
MDI observations as a function of time is expressed by the 
Fourier transform; 
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The same applies to the acoustic egression: 
1 
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The "egression power", 
P(r, t) = \H+{r, t)f 
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is used extensively in holography of acoustic sources and 
absorbers. Equation Q is used to produce egression power 
maps, which show compact positive signatures in the spatial 
and temporal neighbourhoods of localized seismic transient 
emitters. The signature of a localized absorber illuminated 
by ambient acoustic noise is a similarly sharp deficit in egres- 
sion power, appearing as a silhouette against a generally 
positive background when rendered graphically. 

In this study P(r, t) is separately derived from compu- 
tations of H+{r,u}) over 2-4 mHz and 5-7 mHz ranges of 
the spectrum. In practice, there are major diagnostic ad- 
vantages to the 5-7 mHz spectrum, as it avoids the much 
greater quiet-Sun ambient noise at lower frequencies, which 
competes unfavourably with acoustic emission into the pupil 
from the flare. Due to a shorter wavelength, the high fre- 
quency band also provides us with waves that have a finer 
diffraction limit. These advantages come at some expense in 
temporal discrimination, as the egression power signatures 
that result are temporally smeared to a minimum effective 
duration of order 
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This smearing operates in both directions in time, meaning 
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Figure 1. Egression power snapshots of AR10720 on 2005 January 15 taken before, during and after tlie fiare and integrated over a 
5-7 mHz frequency band. The top frames show an MDI visible continuum image of AR10720 (left) at 00:00 UT and a magnetogram 
(right) at 00:28 UT. The second row shows continuum intensity differences 30 seconds before and after the time that appears above the 
respective frames, taken from the GONG observatory at Mauna Loa. The bottom three rows show egression power maps before (row 3), 
during (row 4), and after (bottom row) the flare at 3 mHz (left column) and 6 mHz (right column). The annular pupil for the egression 
computations is drawn in the top left panel. To improve statistics, the original egression power snapshots are smeared by convolution 
with a Gaussian with a 1/e-half-width of 3 Mm. Times are indicated above respective panels, with arrows inserted to indicate the location 
of the acoustic source. Colour scales at right and left of row 3 apply to respective columns in rows 3 — 5. Egression power images and the 
continuum images are normalized to unity at respective mean quiet-Sun values. At 3 mHz this is '~2 kW m~^. At 6 mHz it is 70 W m~^. 



that the acoustic signature of the flare P(r, t), once the com- 
putation is complete, will invariably commence several min- 
utes before the actual onset of the flare and last for several 
minutes afterward even if the actual acoustic disturbance 
was instantaneous (and no seismic signature at all reached 
the pupil until nearly twenty minutes after the flare began). 



4 RESULTS AND ANALYSIS 
4.1 The Seismic Signatures 

AR10720 was predominated by a single (5-configuration 
sunspot. The top row of Fig. [1] shows continuum inten- 



sity (left) and line-of-sight magnetic field (right) of the ac- 
tive region shortly before the flare. The 2005 January 15 
solar flare in AR10720 was classified as XI. 2, localized at 
N14E08 on the solar surface. The GOES satellite measured 
a 1.2 X 10~^Jm~^ X-ray flux in the 1-8 A range integrated 
over the duration of the fiare. Excess X-ray emission be- 
gan at 00:22 UT, reaching a maximum at 00:43 UT, and 
ending at 01:02 UT. There was significant white-light emis- 
sion with a sudden onset, as indicated by the intensity dif- 
ference signatures shown in the second row of Fig. [TJ and 
this coincided closely with hard X-ray (HXR) signatures in- 
dicating high-energy particles accelerated into the chromo- 
sphere. However, unlike the fiares of 2003 October 28-29 
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Figure 2. Acoustic power snapshots of AR10720 on 2005 January 15. Details are the same as for Fig.[T] but local acoustic power maps 
appear in the bottom three rows in place of egression power maps. 



IIDonea fc Lindsevll2005h . there were no sienatures to indi- 
cate the inclusion of high-energy protons in these particle 
influxes. 

To assess seismic emission from the flare, we computed 
the egression, H+, as prescribed by equation (|3} over the 
neighbourhood of the active region at one-minute intervals 
in t, mapping the egression power, P, as prescribed by equa- 
tion (U}, for each value of t. We call a map of P evaluated 
at any single t an egression power "snapshot." From this 
point will refer to the 5-7 mHz bandpass simply as 6 mHz 
and to the 2-4 mHz bandpass as 3 mHz. Egression power 
snapshots before, during and after the flare are shown in 
the bottom three rows of Fig. [1] at 3 mHz (left column) and 
6 mHz (right column). In these computations the pupil was 
an annulus of radial range 15-45 Mm centred on the focus 
(Fig.IBi). 

AU egression power snapshots mapped in Fig. [1] show 
considerably suppressed acoustic emission from the mag- 
netic region, attributed to stro ng acoustic absorption by 
magnetic regions, discovered by iBraun. Duvall fc LaBontg 



lll988t) (see also iBraunI iBraun et al] (|l998h : 

iBraun fc LindsevI |l999a)). Furthermore, all 6 mHz egres- 
sion power snapshots in Fig. [2] also show acoustic emission 
"halos," i.e. significantly enhanced acous tic emission from 
the outskirts of complex active regio ns l|Braun fc LindsevI 
Il999bl : iDonea. Braun fc Lindsevll 19991 ) . 

A conspicuous seismic source is seen in the 6 mHz egres- 
sion power snapshot at 00:42 UT, whose location is indicated 
by an arrow in all of the frames. A close examination of 
the source shows that it has two components. By far the 
most conspicuous component is an intense, compact ker- 
nel ~10 Mm located on the penumbral neutral line of the 
5-configuration sunspot. Somewhat more diffuse but clearly 
signiflcant is a secondary, somewhat lenticular signature dis- 
tributed along the neutral line out to ~15 Mm east and 
~30 Mm west of the kernel. These signatures correspond 
closely with other compact manifestations of the flare. The 
kernel accounts for approximately 45 per cent of the egres- 
sion power integrated over the region encompassing the flare 



The solar flare-induced sunquake of 2005 January 15 5 



soaa 





■^c^u'stJc ro">rr 3 mHi 




D 40 £G E£i TD RO »1 ■ DC- 1 10 1 

TiJlal Iffui^nct 




Kernel Irrmi'ianr'e' 




m ^0 Ml 



Figure 3. Time series of the 3 and 6 mHz egression and acoustic power (integrated over the neighbourhood of the egression power 
signatures) are plotted in the top four rows. The dashed vertical lines mark the time of maximum acoustic emission (00:41 UT) at 
6 mHz. The relatively extended duration of the acoustic signatures is a result of limits to temporal resolution imposed by truncation of 
the spectrum (see equation [5]l . The bottom two rows show visible continuum irradiance at 1 au from the flaring region along the neutral 
line in the neighbourhood of the flare. The emission from the neighbourhood of the kernel component of the 6 mHz acoustic source (plot 
f) is discriminated from the total (plot c). 



signature, with the lenticular component outside of the ker- 
nel accounting for the rest. 

The 3 mHz egression power snapshots shown in the left 
column of Fig. [1] actually show a considerably stronger seis- 
mic emission signature than the 6 mHz signature (right col- 
umn). But, because of the much greater ambient acoustic 
emission at this frequency, the 3 mHz signature is not nearly 
as conspicuous or significant as the 6 mHz signature. It ap- 



pears to have only a diffuse lenticular component and no 
conspicuous kernel to match the 6 mHz kernel. 

It is important to distinguish between the egression 
power, |i?+(r, and the local acoustic power, P(r, t), 

which is the square modulus, \ipiv, t)\^ , of the local wave am- 
plitude i/) at the focus, r. Each pixel in a local acoustic power 
map represents local surface motion as viewed directly from 
above the photosphere. Each pixel in the egression power 
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Table 1. Energy estimates of the seismic signatures of sunquakes detected to date 



Date Class 3 mHz 6 mHz 1 — 8 A X-Rays Visible 

(ergs) (ergs) (ergs) (ergs) 



1996 Jul 09 


X2.6 


7.5 


X 


1027 


2.4 


X 


1026 


2.8 


X 


1029 








2001 Sep 09 


M9.5 


1.1 


X 


1027 


2.0 


X 


1026 


6.2 


X 


1028 


1.2 


X 


1030 


2003 Oct 28 


X17.2 


4.7 


X 


1027 


9.4 


X 


1026 


5.0 


X 


1030 








2003 Oct 29 


XIO.O 


9.4 


X 


1026 


3.5 


X 


1026 


1.5 


X 


1030 


3.8 


X 


1029 


2005 Jan 15 


XI. 2 


2.4 


X 


1027 


1.0 


X 


1027 


3.4 


X 


1029 


2.0 


X 


1030 



map computed by subjacent vantage holography of the sur- 
face is a coherent representation of acoustic waves that have 
emanated downward from the focus, deep beneath the solar 
surface, and re-emerged into a pupil (see diagram of annulus 
in Fig. [TJi) 15-45 Mm from the focus. 

Fig.[2]shows local acoustic power snapshots of AR10720 
at 3 mHz (left column) and 6 mHz (right column) before, 
during, and after the flare. As in the case of egression power 
(Fig. [1} all of the local acoustic power maps show a broad 
acoustic deficit marking the magnetic region. An enhanced 
local acoustic power halo surrounding the active region is 
clearly apparent in the 6 mHz snapshots. The acoustic signa- 
ture of the flare is also clearly visible at 6 mHz. This appears 
to consist of a pair of kernels, a relatively stronger one nearly 
coinciding in location with, but slightly east of, the 6 mHz 
egression power kernel and a weaker one ~10 Mm to west 
and slight north, lying on the neutral line along which the 
lenticular component of the 6 mHz egression power is dis- 
tributed. As in the corresponding egression power snapshot, 
the 3 mHz local acoustic power snapshots show a stronger 
but still less conspicuous signature than that at 6 mHz due 
to a similarly much greater background acoustic power at 
3 Mm. 

Fig. 13] shows plots of the egression and acoustic power 
time series in the 3 and 6 mHz bands and continuum emis- 
sion in the neighbourhood of the seismic signature, discrim- 
inating continuum emission in the region of the kernel com- 
ponent in the 6 mHz egression power signature from the 
total. The flare irradiance profiles were extrapolated by ap- 
plying the assumption that the irradiance is directly propor- 
tional to the GONG continuum signature in the neighbour- 
hood of Ni I 6768 A((d onea fc Lindsevl[2005l ). 

The fiare of 2005 January 15 produced the most con- 
spicuous acoustic signature of any fiare that has produced 
a detectable seismic emission. This appears to be because 
such a large fraction of the energy was released into the 
high-frequency (5-7 mHz) spectrum, where the competing 
ambient acoustic power is so far suppressed. Table |4T1 shows 
the energy estimates of the seismic transients radiated into 
the active region subphotosphere by five fiares that have pro- 
duced conspicuous seismic signature^ compared with en- 

^ The energy estimates in Table l4!T1 were obtained by integrating 
the egression power over the neighbourhood of the seismic sources 
(e.g. those shown in Figs. lg,h for the 2005 January 15 flare). This 
computation is blind to waves that miss the 15-45 Mm in the 
first skip. Comparative seismic holography applied to simulated 



ergy emitted in X-rays in the first 20 minutes of the fiare. 
It should be noted that the 3 mHz energy for the flares 
preceding the 2005 January 15 flare are actually calculated 
at 3.5 mHz. Extrapolating through the missing 4-5 mHz 
acoustic spectrum for the flare of January 15, we project a 
total acoustic emission of ~ 4 X 1020 J (~ 4 X 1Q27 erg). 

4.2 Visible Continuum Emission 

Various aspects of visible continuum emission during the 
2005 January 15 flare are shown in Figs. [T] [2] and [S] The 
visible-continuum images in Figs. [T] and [5] were obtained by 
MDI at 00:00 UT, ~37 min before the onset of the flare. 
We obtained visible continuum maps of AR10720 during 
the flare from the GONG observatory at Mauna Loa. Tech- 
nically, the GONG "continuum intensity maps" represent a 
measure of radiation in a ~1 A bandpass centred on the 
Ni I 6768 A line, whose equivalent width is only a fraction 
of an A . Frames c) and d) in Figs. [1] and [2] show the dif- 
ference in continuum intensity between the GONG images 
30 seconds before and after at the time indicated above the 
frame. Continuum emission is elongated along the magnetic 
neutral line, corresponding closely to the lenticular compo- 
nent of seismic emission seen at 00:42 UT in Fig. [T]i. The 
brightest emission seen in the intensity difference shown in 
Fig. [ijl comes from a very compact kernel whose location 
coincides very closely with that of the conspicuous kernel of 
6 mHz emission (Fig. [l]i). 

If we assume that the continuum emission emanates 
isotropically from an opaque surfac^El the resulting esti- 
mate of the total energy emitted in the visible continuum is 
2.0 X 1023 J (2.0 X 1030 erg). This is ~500 times the total 
seismic energy we estimate the fiare to have emitted into 

acoustic transients, and to MDI observations of flares with dif- 
ferent sized pupils, indicate that the energies quoted in Tabic IXTI 
account for 80 — 95 per cent of the total, depending on the source 
distribution. 

2 The assumption is that the specific intensity is independent of 
the vantage, which implies that the total intensity decreases in 
proportion to the cosine of the angle of the vantage from normal 
as a result of foreshortening. If the source was assumed to be 
optically thin, the resulting energy estimate would be greater by a 
factor of two. The former appears to be the more realistic estimate 
for the fraction of visible continuum radiation coming directly 
from the chromosphere, based on the thesis that ionization of 
chromospheric hydrogen at the temperature minimum renders the 
low chromosphere opaque. 
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Figure 4. MDI Doppler-difTerence images showing the expanding ring-shaped wave packet produced by the 2005 January 15 flare. The 
times shown above respective panels represent one-minute difi'erences. The arrows pointing in the South-East direction (i.e. upper arrows) 
show the location of the upper-half of the wave front while arrows pointing in the North- West direction (i.e. lower arrows) indicate the 
lower half of the wave front. Grey-scale at top left expresses Doppler velocity differences in units of ms~^ and applies to all frames in 
the figure. 



the holographic pupil. Continuum radiation into the neigh- 
bourhood of the 6 mHz kernel signature was 6.0 x 10^^ J 
(6.0 X 10^^ erg). This accounted for ~30 per cent of the 
total, as compared to 45 per cent of the 6 mHz seismic sig- 
nature. Continuum emission from in the neighbourhood of 
the 6 mHz kernel was significantly more sudden than that 
of the remainder of the acoustic signature. 



ings that relatively small fla res can emit dispr oportionate 
amounts of acoustic energy jPonea et aLlliooil ). However, 
even in these cases the fraction of the energy that is re- 
leased by the flare into the solar interior acoustic spectrum 
remains relatively small. 



The 2005 January 15 flare contributes to recent find- 
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Figure 5. MDI Doppler-difTerence images showing ring-shaped wave packet accelerating outward from the site of the 2005 January 
15 flare. The left panel shows the local Doppler signature along the magnetic neutral line during the impulsive phase of the flare at 
00:40:30 UT. Arrows in the right panel indicate locations of the ripples propagating outwards from the site of local disturbance 25 
minutes later. 



4.3 The Surface Ripples 

Holography allows us to image the acoustic source of the 
sunquake when the surface manifestation of the seismic 
emission is difficult to detect. In the case of the exception- 
ally powerful seismic transient from the flare of 2005 Jan- 
uary 15, the surface signature is quite evident in the raw 
MDI Doppler observations, a point to which A. Kosovichev 
(2005, private communication) drew our attention shortly 
after we reported the discovery of the sunquake to him. To 
extract the seismic oscillations in the observations we sub- 
tracted consecutive MDI Doppler images separated by one 
minute in time. We applied this Doppler-difference method 
to a period of observation (~1 hour) around the time of the 
flare. Results are shown in Figs. |4] and (5] 

The Doppler signature of the flare is clearly evident at 
00:40 UT (Fig. O left panel). At approximately 20 minutes 
after the appearance of the flare signature in the sunspot 
photosphere (at 01:00 UT), we are able to see the seismic 
response of the photosphere to the energy deposited by the 
flare in the form of "ripples" on the solar surface. In the 
sequence of one-minute Doppler-difference images in Fig. |4j 
we can see the asymmetrical ring-shaped wave packet prop- 
agating from the flare site with the first wave-crest appear- 
ing approximately 12-15 Mm from the flare in a North- 
Easterly direction. The lower half of the wave-packet has 
a much smaller amplitude and is propagating in a South- 
Westerly direction. The arrows in Fig.|3]indicate the location 
of the observed wave fronts. The Doppler-difference images 
in Fig. [5] show a close-up of the active region at the time 
of the flare (at 00:40:30 UT, left panel) and the resulting 
ring-shaped wave packet (at 01:05:00 UT, right panel). 

The wave-packet was seen to propagate to a maximum 
distance of approximately 21 Mm from the flare signal, hence 
travelling a total distance of 6 - 9 Mm and lasting for about 
8 minutes on the surface, after which the wave amplitude 
dropped rapidly and the disturbance became submerged in 
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Figure 6. TRACE white-light image for AR10720, on 2005 Jan- 
uary 15 (00:17:54 UT) with the 12-25 keV RHESSI contours (10, 
20, 30, 40, 50, 60 and 80 per cent of the maximum flux). The (0,0) 
coordinates correspond to the location of the seismic source. 



the ambient noise. The lower half of the wave-packet (propa- 
gating towards the South- Western part of the active region, 
indicated by the lower of the two arrows superimposed on 
the Doppler-difference images in Figs. |4] and [5]) was much 
smaller in amplitude and obscured for much of the 8 min- 
utes. 
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Figure 7. The HXR flux in tlie 4-25 (0.5-4 A) and 25-50 keV energy ranges as observed by GOES (dot-daslied curve, multiplied by 
a factor of 10®'^ Wm~^) and RHESSI (bold curve, multiplied by a factor of 10~^'* counts) respectively. The solid curve represents the 
1-minute mean averages of the Doppler velocity (ms~^) in the quake region for the period 00:00:00-01:20:00 UT. The dashed vertical 
line represents the observed maximum emission at 6 mHz. 



4.4 Hard X-Ray Emission 

The TRACE data for the 2005 January 15 flare in the 
white-light channel have a variable cadence for the period 
00:00:00-01:00:00 UT. Fig.|6]shows the TRACE white-light 
image taken at 00:17:54 UT, approximately 10 minutes be- 
fore the onset of the XI. 2 flare with the RHESSI 12-25 keV 
contours overlaid. The RHESSI HXR image is averaged over 
the period 00:41:33-00:42:34 UT. The time of peak intensity 
in this energy band occurs at 00:42:04 UT, a close temporal 
correlation with the maximum of the seismic emission de- 
tected at 6 mHz. The HXR emission is thought to represent 
bremsstrahlung emission from hig h-energy coro nal electrons 
impinging into the chromosphere (|Brownlll97ll ). 

The 12-25 keV emission at 00:42:00 UT extends along 
the neutral magnetic line. We identify three compact HXR 
sources (see the numbers in Fig.[S|. Source 2 is the strongest, 
while source 3 is the weakest. These could represent the 
foot-points of a complex magnetic loop. However, source 1 
(which emits 50 per cent of the total flux) spatially coincides 
with the lenticular component of the 6 mHz seismic source 
(see Fig. [1]). This reinforces the role of non-thermal parti- 
cles in supplying the energy that drives the seismic emis- 
sion. Similar comparison s have been observed in other flares 
l|Donea fc Lindsevll2005l ). 



Furthermore, Fig. [7] reveals that the velocity impulse of 
the flare in the sunspot photosphere was almost as sharp as 
the HXR flux detected in the 4-25 keV (0.5-4 A) energy 
range by the GOES satellite, but the maximum HXR emis- 
sion (observed at ~00:43:00 UT) appears to have occurred 
~2 minutes after the maximum velocity depression at the 
photosphere (00:41:00 UT). In fact, a sudden drop of approx- 
imately 100 ms~^ in the mean velocity of the Doppler signal 
(an upflow) is observed in the 3 minute period from 00:38:00- 
00:41:00 UT. The RHESSI HXR peak in the higher energy 
band of 25 - 50 keV plotted in Fig. |7]occurs at ~00:41:00 UT, 
which temporally coincides with both the maximum of the 
seismic source at 6 mHz and the velocity depression at the 
photosphere. We also note that the peak emission in the 
3-12 keV energy band detected by both GOES (1-8 A) 
and RHESSI (occurring at ~00:44:00 and ~00:47:00 UT 
respectively, but not plotted), also have a close temporal 
correlation with the maximum of the seismic emission. 



5 DISCUSSION AND SUMMARY 

The X1.2-class flare of 2005 January 15 produced one of 
the most powerful sunquakes detected to date and by far 
the most conspicuous, on account of exceptionally powerful 
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emission above 5 mHz from a compact source. Certain qual- 
ities exhibited by the flare of January 15 are shared by all 
other known acoustically active flares. The first is the coinci- 
dence between strong compact acoustic sources and nearby 
signatures of HXR emission. This suggests that high-energy 
particles supply the energy that drives the acoustic emis- 
sion, and it is evident from the electromagnetic emission 
attributed to these particles that they contain more than 
sufBcient energy for this purpose. The appearance of sud- 
den, conspicuous white-light emission from the flare of 2005 
January 15 closely co-spatial with the location and morphol- 
ogy of the holographic signatures is similarly characteristic 
of a h other known acoustical l y acti ve flares so far. 

iKosovichev fc Zharkoval (|l998h proposed that seismic 
emission into the solar interior in sunquakes is the con- 
tinuation of a chromospheric shock and condensation re- 
sulting from explosive ablation of the chromosphere and 
propagating downward through the photosphere into the 
underlying solar interior. Chromospheric shocks are well 
known under such circumstances, based on red-shifted Hq 
emission at the flare site at the onset of the flare. The 
theory of their dynamics was worked out at length by 
iFisher. Canfield fc McClvmond l| l985allblH) and others since. 
The hypothesis that solar interior emissio n is a direct contin- 
uation of such shocks was considered by iDonea fc LindsevI 
l|2005l l, who found the signature of a strong downward- 
propagating chromospheric transient in Na Di-line observa- 
tions of the flare of 2003 October 29. However, we are now 
aware of similar chromospheric transient signatures with no 
significant attendant holographic signature to indicate seis- 
mic emission into the active region subphotospher^fl 

In these instances, the signatur e of sudden white-ligh t 
emission is relative ly weak. Following iMachado et al.l(|l989t ). 
iDonea et al.l (|2006l ^ proposed to attribute the lack of seismic 
emission where there is a strong chromospheric transient 
but only a weak or absent white-light signature to strong 
radiative damping that depletes the chromospheric transient 
before its arrival into the low photosphere. 

In all acoustically active flares encountered to date, 
there is a strong spatial correlation between the sources 
of seismic emission and sudden white-light emission. This 
remains conspicuously the case for the flare of 2005 Jan- 
uary 15, as a comparison between Figs. Id and Ih shows. 
In some instances, e.g. the large flares of 2003 October 29, 
the source of the white-light emission has been much more 
extensive than the source of the acoustic emission, the for- 
mer many times the area of the latter and encompassing it. 
However, in these instances the temporal profile of visible 
continuum emission significantly away from any of the sites 
of seismic emission has been comparatively sluggish and dif- 
fuse. What has particularly and consistently distinguished 
the white-light signature in the neighbourhood of the acous- 
tic emission has been the suddenness of its appearance, on 



^ An example is seen in the western foot-point of the magnetic 
loop that hosted 2003 October 29 flar e (see right frame in secon d 
row of fig. 8 labelled "Red [0 min]" in lDonea fc LindsevI ||2005|) ). 
The corresponding signature of sudden white-light emission, seen 
in fig. 9 of the same, shows only a weak signature at the same 
location. The seismic signature, seen in the lower left corner of the 
upper right frame of fig. 11 of the same, shows correspondingly 
weak seismic emission. 



a time scale of a minute or two, and possibly considerably 
less than a minute given that the observations of continuum 
emission associated with flares to date have been limited to 
a cadence of one minute. 

It should be kept in mind that the energies released 
in known seismic transients have invariably been a small 
fraction of the energy released into the visible continuum 
spectrum. The actual fraction has varied considerably, from 
a few millionths, in the case of the flare of 2003 October 29 
(|Donea fc Lindsevll2005l ) . to a few thousandths, for the flare 
of 2005 January 15. However, if only the sudden-onset con- 
tinuum emission in the neighbourhood of the seismic sources 
is included, then the ratio for the flare of 2003 October 29 is 
similar to that of the flare of 2005 January 15. This is what 
is listed in Table |4T1 of this study. 

The close coincidence between the locations of sudden 
white-light emission and seismic transient emission in all 
acoustically active flares to date suggests that a substantial 
component of the seismic emission seen is a result of sud- 
den heating of the low photosphere associated with the vis- 
ible continuum emission seen. A complete analysis of wave 
emission as a result of transient heating involves detailed 
considerations of energy and momentum bal ance. An ap- 
proxim ate account of these was undertaken bv lDonea et al.l 
l|2006t ). Basic considerati ons of mome n tum b alance are de- 
scribed in Section 4.3 of 'Ponca e t al.l l|2006l ), adapting the 



discussion by ICanfield~et al.. (,1990l ) of momentum balance 
in chromospheric transients to transients similarly excited 
by s udden heating in th e low photosphere. 

IDonea etH1 120061) devised a rough, preliminary phys- 
ical model to estimate the energy of the seismic transient 
to be emitted as a result of sudden, momentary heating of 
the low photosphere to a degree consistent with the tran- 
sient white-light signature closely coincident with the seis- 
mic source in the M9.5-class flare of 2001 September 9. Their 
estimate expressed the energy, E, of the seismic transient in 
terms of the thermal energy, U, radiated or dissipated into 
the low photosphere, and the fractional increment, Sp/p, in 
pressure that would result from the heating: 



E 



2 p 



1 5p 

3 7 
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(6) 



where H is the e-folding height of the photospheric density. 
This relation appears to be roughly consistent with the few- 
percent continuum intensity variations observed for the flare 
of 2001 September 9, if the relation between 5p and 51, the 
variation in continuum intensity, can be approximated by 
the Stefan-Boltzmann law. 



5p 
P 



ST 



15/ 
iT' 



(7) 



and the heating is accomplished within a duration not ex- 
cessively longer than Tac ~ l/i^ac ~ 40 s, where uiac is the 
aco ustic cutoff frequency in t he low photosphere (see Section 
4 of lDonea fc LindsevI l|2005h ). A similar exercise applied to 
the flare of 2005 January 15 leads to similar results. In fact, 
the ratio of the seismic energy to the electromagnetic en- 
ergy is roughly the same for both of these flares, as are the 
mean intensity increments if credible boundaries are chosen 
over which to take the mean. To the extent that we can re- 
solve the fine details, acoustic emission from the flare of 2005 
January 15 could reasonably be the result of photospheric 
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heating similar to that of the 2001 September 9 flare but 
over approximately twice the area. Differences between the 
two flares could be attributed to differing photospheric or 
subphotospheric thermal conditions and differing magnetic 
fields, for which the foregoing approximation contains no 
account. 

A detailed examination of the physics of heated mag- 
netic photospheres is needed to lend credibility to the hy- 
pothesis that seismic emission from acoustically active flares 
is driven by sudden heating of the low photosphere by any 
mechanism whatever. At this point we will only say that 
this hypothesis appears to be consistent with our present 
limited understanding of the observations. However, there 
is some controversy as to the implications of visible contin- 
uum emission during flares with respect to heating of the 
low photosphere. In the case of the flares of 2003 October 
28 - 29, the signature of high-energy protons along with the 
particles that gave rise to X-ray emission lent considerable 
weight to the interpretation of visible continuum emission 
in terms of a heated low photosphere, as protons are suf- 
ficiently massive to penetrate to the bottom of the photo- 
sphere and heat it directly by collisions. The flare of 2005 
January 15, on the other hand, confronts us with an instance 
of intense seismic emission with no indication of high-energy 
protons among the energetic particles that supply the en- 
ergy on which the acoustic emission depends. Energetic elec- 
trons consistent with HXR signatures cannot penetrate into 
the low photosphere in anywhere near sufficient numbers to 
account for the heating required by the seismic signa tures 
ijMetcalf. Canfield fc Sa ba 19901. IChen fc Dind (|2006l ) also 
affirm that the white-light fiare signatures highlight the im- 
portance of radiative back-warming in transporting the en- 
ergy to the low photosphere when direct heating by beam 
electrons is impossible. 

In such cases, it appears to be well established 
that the origin of white-light emission would have to 
be entirely in the chrom osphere, where energetic elec- 
trons dissipate their energv (iMetcalf. Canfield fc Saballl990l : 
IZharkova fc Kobvl inskii"l99ll. ll993l ')~ mainly by ionizing pre- 
viously neutral chromospheric hydrogen approximately to 
the depth of the temperature minimum. Nevertheless, even 
in these instances, it appears that the low photosphere itself 
would be significantly heated as a secondary, but more or 
less immediate, effect of chromospheric ionization. This is 
primarily the result of Balmer and Paschen continuum edge 
recombination radiation from the overlying ionized chromo- 
spheric medium, approximately half of which we assume 
radiates downward and into the underlying photosphere. 
When the intensity, 51, of this downward flux is commensu- 
rate with a temperature perturbation, 5T, consistent with 
the Stefan-Boltzmann law (equation [7] ) , the result of such a 
flux is understood to be heating of the low photosphere such 
as to bring about a tempe rature increment of roughly this or- 
der within a few seconds (jPonea et al ] |2006l : lMachado" et al.l 
1 19891 : iMetcaff et "all 120031 ) . Heating of the photosphere by 
the mechanism desc ribed above is known as "back-warming" 
l|Metcaff et al.ll2003h and a substantial fraction of the contin- 
uum emission seen in white-light flares is thought to repre- 
sent the downward flux from an ionized chromosphere ther- 
mally re-emitted by the now heated photosphere. In this 
light, the strong correlation between sources of white- light 
and seismic emission into the solar interior might be re- 



garded as strong support for the back-warming hypothesis 
when this relation persists in flares devoid of protons among 
the high-energy particles that drive the flare. This is cer- 
tain ly the case for the flare of 200 5 January 15 . 

iDonea fc LindsevI l|2005h and lDonea et al.1 (|2006l ) sum- 
marize our understanding of the relationship between the ef- 
ficiency of seismic emission and the suddenness of the heat- 
ing that drives the seismic transient. Based on these con- 
siderations, one has to suspect that the perceptibly more 
sudden profile of continuum emission in the neighbourhood 
of the kernel component of the 6 mHz emission accounts to 
a significant degree for the disproportionate power in the 
6 mHz egression-power signature. This is one of the many 
aspects of fiare acoustics that would benefit from detailed 
modelling, including a careful account of magnetic forces. 

At this point, our understanding of seismic emission 
from flares remains relatively superficial. However, evidence 
for the general involvement of photospheric heating is now 
considerable. What is needed for further understanding is 
detailed modelling with a careful account of the physics, 
including radiative transfer and magnetic forces in realis- 
tic sunspot photospheres and subphotospheres. With such 
an understanding, acoustic emission from flares could con- 
tribute major benefits to seismic diagnostics of active region 
subp hotospheres and the physics of mode conversion (|Callvl 
I2OOOI) . 
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